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AN  ERROR  ANALYSIS  FOR  THE  FINITE  ELEMENT  METHOD 


APPLIED  TO  CONVECTION  DIFFUSION  PROBLEMS 


ABSTRACT 

This  paper  analyzes  the  finite  element  method1  applied  to  a  convection 

diffusion  model  problem.  Linear  elements  are  used  for  the  trial  space.  The 

error  is  measured  in  a  norm  closely  related  to  the  L  norm.  When  the  test 

P 

space  is  composed  of  linear  elements  with  parabolic  upwinding,  the  method  is 
shown  to  be  optimal  when  the  input  data  is  piecewise  smooth  —  a  condition  which 
is  usually  observed  in  practice.  Without  these  smoothness  assumptions,  the 
method  is  shown  to  be  non-optimal,  even  if  the  class  of  test  spaces  is  extended 
to  include  any  elements  which  have  a  shape  independent  of  the  mesh  size. 

\ 

\ 


INTRODUCTION 


A  large  amount  of  attention  has  recently  been  focused  on  the  application 
of  finite  element  methods  to  singularly  perturbed  boundary  value  problems. 

These  problems  arise,  for  example,  in  convection  diffusion  equations  for  fluid 
mechanics  in  which  the  convective  term  dominates.  Finite  element  methods 
which  employ  the  standard  piecewise  polynomial  test  and  trial  spaces  lead  to 
solutions  having  spurious  oscillations  unless  the  mesh  size  h  is  excessively 
small.  A  popular  way  to  alleviate  this  problem  is  to  "upwind"  the  test  space. 

In  the  case  of  linear  elements  this  can  be  done  by  adding  a  quadratic  term, 
multiplied  by  some  parameter  n  «  to  each  piecewise  linear  basis  function  of 
the  test  space  (see  e.g.  [5],  [10]-[12],  [ 15 ] —[17]).  Heinrich, et  al.  [11],  [12], 
Christie,  et  al.  [6]  and  others  have  displayed  the  optimal  a  ,  which  is  deter¬ 
mined  in  such  a  way  so  that  for  the  problem 

-cu"+u'  =  0  ,  u(G)  =  0  ,  u ( 1)  =  1  , 

the  approximate  solution  agrees  with  the  exact  solution  at  the  nodes.  Griffiths 
[10]  selects  the  parameter  a  in  order  to  obtain  the  best  quasi-optimal  type 
estimate  of  the  form 


I  I u-uh  1  I  -  c<e»h)  I  !u-w|  | 

w€S 

n 

Where  u  is  the  exact  solution,  u^  is  the  approximate  solution,  is  the 

trial  space,  and  ||*||  is  the  norm  of  the  derivative.  However,  the 

constant  c  is  not  bounded  uniformly  in  e  and  h  ;  for  example, if  c  is 
very  small  with  respect  to  h,  c(e,h)  will  increase  with  rate  h'1  as  h  -  0  . 

Upwinding  of  the  test  space  can  also  be  done  through  the  use  of  L-spline 
basis  functions.  Methods  using  these  spaces  are  studied  in  Hemker  [13], 
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de  Groen  [7],  and  Hemker  and  de  Groen  [8],  but  a  norm  is  used  to  measure  the  error 
between  the  exact  solution  u  and  its  piecewise  linear  approximation  in  which 
the  error  cannot  be  small  unless  h  <<  e. 

The  analysis  in  this  paper  is  done  specifically  for  the  model  problem 

(1.1)  -eu"  +  u’  =  f  , 

u(0)  =  u(l)  =  0 

In  section  2  we  present  the  general  theory  for  obtaining  quasi-optimal  estimates,  as 
established  in  Babuska,  Aziz  [1].  Mesh  dependent  norms  for  the  model  problem, 
analogous  to  those  presented  in  [2],  are  defined  in  section  3.  We  will  be  measur¬ 
ing  the  error  in  a  norm  similar  to  the  norm, which  is  appropriate  in  our 

case  because  it  allows  for  approximability  when  using  a  piecewise  linear  trial 
space.  In  section  4  a  quasi-optimal  result  of  the  form 

(1.2)  l|u-u.||_<c  inf  )  |  u— w  j  |  , 

w€Sh 

where  c  is  independent  of  e  and  h  ,  is  proven  when  the  test  space  is  composed 
of  L-spline  basis  functions.  However,  it  is  also  shown  that  for  a  large  class  of 
test  functions,  it  is  in  general  impossible  to  obtain  a  result  of  the  form  (1.2). 

In  particular,  quasi-optimality  cannot  be  proven  with  the  use  of  quadrati cally 
upwinded  elements  for  any  choice  of  the  parameter  a  unless  we  restrict  our¬ 
selves  to  the  case  of  "reasonable"  inputs  f  .  This  is  shown  in  our  main  result, 
presented  in  section  5,  which  says  that  if  f  is  piecewise  smooth,  and  a  is 
chosen  to  be  the  optimal  one  introduced  by  e.g.  Christie  et  al.  16],  and 
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Zienkiewicz  and  Heinrich  [17],  then  quasi-optimality  is  obtained. 


2.  SOME  ABSTRACT  RESULTS 

In  this  section  we  review  two  crucial  results  concerning  variationally 
formulated  boundary  value  problems  and  finite  element  approximations. 

Theorem  2.1.  Let  ^  and  ^  A  be  two  reflexive  Banach  spaces,  indexed 
by  a  parameter  A  with  A  varying  over  some  index  set,  with  norms  j  |  • j  |  _  and 
|  |-  |  1 2  A  respectively,  and  let  B  be  a  bilinear  form  on  &  *  K2  A  "  We 
suppose  the  following  are  satisfied: 


(2.1) 


Iba(u’v> I  i  cxl  M 


for  all  u  ^  , 


v€K 


2, A 


(2.2) 


inf 

u€K 


1,A 

!lfA 


=1 


sup  |B.(u,v)  |  >_  C7  >  0 

V€K2,A 

M  l2>A*l 


and 


(2.3) 


sup  |b.(u,v))  >  0  ,  for  each  0  4  v€K, 

u€K  . 

1,A 


where  and  C2  are  positive  constants  ,  possibly  depending  on  A  . 

Then  if  f  6  (^  ^)  ’  ,  there  exists  a  unique  solution  u€K^  ^  to  the  problem 


(2.4) 


Ba(u,v)  -  f(v)  ^v€K2,A 
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Moreover,  u  satisfies 


M  1 3., a  —  c2  I  lf  I  Ik; 


2,  a 


If  the  bilinear  form  B  (•»•)  satisfies  the  assumptions  (2.1),  (2.2),  and 
(2.3),  B^  is  said  to  be  a  (C^,C2)-proper  bilinear  .orm  over  the  space 
Kl  ^  x  K,  j  .  It  should  be  noted  that  (2.2)  and  (2.3)  can  be  shown  to  be 
equivalent  to 


(2 . 2) : 


inf 

veK2.4 

|v||2,a-1 


sup  |b  (u,v)  |  Cy  >  0 


u€  K 


1,A 

1,A=1 


and 


(2.3) 


v€K 


sup  |Ba(u,v)(  >  0  ,  Vo  +  uC^  t 


2, A 


This  observation  will  be  specifically  used  in  this  paper. 

Since  we  will  be  studying  finite  element  approximations  to  u  ,  we  let  A 

and  S2  A  be  finite  dimensional  subspaces  of  A  and  ^  A,  respectively. 

Clearly  condition  (2 . 1)  holds  on  S.  .  *  S,  .  with  the  same  constant  C.  .  We 

1,A  2, A  1 

will  be  invoking  the  following  theorem  concerning  the  finite  element  solution  u^ 


Theorem 

2.2. 

Suppose  is 

(C^ , cp-proper  over 

.  *  .  furnished 

A  2,  A 

with  the  norms  | 

N llfA  and  |  |- |  1 

2  respectively.  Let 

u^K,  ,  ,  and  let 

1  *  A 

Uh€Sl,A  be 

the 

unique  solution  to 

BA(uh,v)  =  B&(u,v)  for 

all  v€S2  a  •  Thi 

The' 
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3.  MESH  DEPENDENT  NORMS  AND  SPACES 

In  this  section  we  define  the  various  norms,  spaces,  and  bilinear  forms  used 
throughout  this  paper.  The  norms  introduced  here  are  analogous  to  those 
defined  in  [2]. 


Throughout  the  paper  Hp(I)  ,  k  =  0,1,...,  1  <  p  <  »  ,  will  denote  the  usual 


Sobolev  space  on  the  interval  I  in  R1  consisting  of  functions  with  k 


derivatives  in  L^(I)  .  On  this  space  we  have  the  usual  norms  given  by 

f  -  r  f  i  ( j )  , 


(3.1) 


'k,p,I 


(  \  |u(Jdx)  |Pdx)P  ,  1  <_  p  < 


j=0 


J  ess  sup  |u(j) |  »  p  = 

j=0 


•1  1 

Hp (I)  denotes  the  subspace  of  H^(I)  of 

of  I  .  Note  that  H°  *  L 

P  P 

Let  A  =  (0  =  xq  <  x1  <  ...Xjj  *  1}  , 


on  the  Interval  I  =  [0,1]. 

Let  hj 

=  XJ 

Pj  =  (hj  +  hj+1)/2  ,  j  =  1,. 

.  ..N-l  , 

and 

We  now  define  the  space 

H°  A  ’ 

P ,  A 

1  < 

furnished  with  the  norm 


functions  that  vanish  at  the  endpoints 

where  N  =  N(A)  ,  be  an  arbitrary  mesh 

'  xj_i  ’  Xj  =  (xj-l,xj)  •  J  = 

h(A)  =  max  h.  . 

j  J 

o  1 

p  ,  to  be  the  completion  of  H  (I) 


(3.2) 


IN 


'  1 
( 

0 


p,A 


N-l 


1 

u  |  Pdx  +  J  p  .  |  u(x . )  |  P)P  .  1<^P<”  , 

1  =  1  J  J 


V 


p  = 


The  space  H°  .  can  be  easily  identified  with  I.  ©  1R1 

p,A  J  p v 


,  that  is. 
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u  =  (u.d^  . . .  ,dN_1)  €H°^ a  =  Lp©mN  1  ,  and 


(3.3) 


l“N[  (I) 

P 


A 

+  J  p  .  ]  d  . 
j=l  J  J 


'  i  P  , 

}  ,  1  £  P  <  “  , 


max0|u||L  (I),  |dj  |  }  , 


In  consistency  with  our  definition,  we  say  u€h°  aHH^(I)  if  u£H^(I)  and 

p,A  p  p 

dj  =  n(Xj)  ,  j  =  1 , . . . N-l  . 

Let  us  now  define  ^  =  {v^H^CI):  vj^  €H^(Ij),  j  =  .  We  will 

equip  this  space  with  a  norm  to  be  defined  later. 

On  H°  .  x  .  where  —  +  —  =  1  ,  we  define  a  bilinear  form  B.(-,0 

p,A  q,A  ’  P  q  A 


(3.4) 


IN  f  W-l 

i,v)  =  l  u(-ev"-v’)dx  -  l  ed.J(v’(x.)) 

j  =  ljT  J-l  J  3 


where  J(v'(x,))  =  v'(x.+0)  -  v'(x.-O)  and  v'(x.+0)  =  lim  v'(x)  .  These  limits 
J  J  J  3  x+x^ 

2  ^ 

are  well  defined  because  v[T  €H  (I.)  for  each  k  .  Now  we  will  furnish  the 

l.  q  K 

2  k 

space  11^  ^  with  the  norm  |  |  |  •  |  |  j  ,  defined  by 

|B.(u,v) | 

i  i  i  l  i  l  A 


It  is  evident  that  the  triangle  inequality  and  linearity  are  satisfied  by 

III*  || |  .  We  must  show  positive  definiteness.  To  do  that  we  make  use  of  the  following 


identity  for  v€H 
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(3.6)  v(x) 


X  _  X-t 

(L*v)(t)(l-e  E  )dt  - 

. 

0 


I  (e 


xiix 


(x-x1) 

E  -l)eJ(v’(x.)) 


+ 


(x  -1) 

[n-i  — 

I  (e  e  -l)eJ(v’(x.)) 


1  _  1— x 

(L*v) (t) (1-e  £  ) 

0 


where 


(3.7)  (L*v) lT  =  (ev"+v' ) j T  ,  j  - 

j  j 


For  v€H  A  ,  select  uq  =  ,dj  ,d2 , . .  .dN_j )  such  that  uL  =  sgn(L*v)j 

j  j 

and  d^=  -  sgnJ(v’(x.))  ,  j-l,...,N-l  ,  then 


(3.8) 


N-l 


|B(u  v)|  =  l  | L*v | dx  +  J  e|J(v'(x  )) j  >  0  . 

j=l\  j=l  J 


If  B(u  ,v)  =  0  then  L*v |  =  0  and  -J(v'(x.))  =  0  ,  and  (3.6)  yields 

°  Ak  J 


v  *  0  .  Therefore, 


is  positive  definite,  and  hence,  a  norm. 


Let  us  now  introduce  another  norm.  For  any  v£H  .  define 

p,A 


by 


(3.9) 


£  |ev"+v'|pdx  +  £  ep | J(v' (x, )) |Pp*  p 

j  =  IJ  j-i  J  J  . 


H 


P,t,A 


if  1  £  p  < 


max  | lev 
i<i  <N 


'+v' | j  +  max  e|J(v'(x.)|o 
l<j<N-l  J 

if  P5 


-1 


p,e, 
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We  shall  now  prove  that  the  norms 


for  1  £  p  <  “>  and  equivalent  when  p  = 


Lemma  3.1.  Let  v£H  ,  then 

p,A 


are  equal 


P,£  .A 


(3.10) 


1  <  p  < 


p,e,A 


” ,  r,  A 


”»c,  A 


Proof.  For  1  <  p  <  a  straightforward  application  of  Holder's  inequality 


yields 


f1  -  N  - 

|  B.(u,  v)  |  <_  (  |u|qdx)q(  J  |ev"+v'  |Pdx)P 

i=lj. 


+  (l  p  |d.|q)q(  [  cp|J(v' (x  )) |pp!  P)P 
j=l  3  3  j=l  J  J 


Tr1  N-i  -i-r  n  c 

lJ 


I  cv"+v'  |Pdx 


YrP|J(v'(x  ))|PpJ'PlP  = 
3=1  3  3  J 


Ml  .  M vl  | 

H  .  H  . 
q,A  p,E,A 


and  thus 


p,f  , A 


(3.11) 
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That  3.11  holds  when  p  =  1  and  p  =  00  follows  directly  from  Holder's  in¬ 
equality  . 

2 

For  the  inequality  in  the  other  direction,  let  v€H  ,  1  p  <  “,  be 

p  ,  A 

_  Q 

given.  Select  u  =  (u,d,,...dv,  ,)€H  ,  as  follows: 

V  1  N-l  q ,  A 

u  =  -|l*v|P  ^sgn(L*v)  , 
dj  =  -eP~1|J(v' (Xj)) jP  1oj  Psgn(J(v' (xj)) 

for  1  <  j  <  N-l  . 


Then  u 


V  H°  A 

q,A 


,  1  <  p  < 


H 


p,c,A 


u  =  1  , 

V  H° 

H»  .A 


p=l  , 


and  Ba(uv,v)  =  j  |v| |^2 

H  . 
p,c,A 


,  IIP 


<  OO 


This  immediately  yields 


M 1 1  >  |  M  |  2  ,  i  <  P  <  co 

H  . 
p,c  ,A 


For  p  =  c°  ,  let  L  be  the  index  such  that 


|  t:v"+v'  |  |  (  =  max  ||ev"+v'||  , 

L»UI/  l<j<N  L-UJ) 


and  let  J  be  the  index  such  that 


t  |j(v' (x.))  Ip.1  •=  max  e |J(v' (x.)) |p  1 
1<J<N-1  J  ' 


1] 


Let  n  >  0  be  given  and  define 


En  =  (x€Il:  |(cv"+v')(x)  |  >  |  |ev"+v' |  |L  (J  }  -  n) 


Then  m(E  )  >0.  where  ra(A)  is  the  Lebesgue  measure  of  A. 

n  ’ 

For  v€H^  ^  define  u^  =  (u,d^, . .  .d^_^)  as  follows: 


(3.12) 


u  =  -x„  (m(E  ))  1sgn(L*v)  , 

En  n 


d.  = 
J 


r pJ1sgn(J(v' (Xj)) ,  if  j  =  J  , 


lo 


j  =  1,...J-1,J+1,...N-1  , 


where  x»  denotes  the  characteristic  function  of  the  set  A  .  Then 
A 

I  I u  | I  =2  and 

II  V  1  1  o 

H1  A 

1,A 

I | | v | | j  >  (| |ev"+v' | |L  (  j-n  +  e|J(v'(Xj)) | P / 2  . 

00  AL 


Since  n  was  arbitrary. 


Ml 


H 


1/ 


v  > 


°°,e  iA 


which  concludes  the  proof. 


— ^  That  this  estimate 


e  is  sharp  can  be  seen  by  taking, for  example,  A  =  {O.y.l}  , 
x  .  0  <  x  <  1. . 


c  =  ^  .  v(x)  = 


r  x  .  u  <  x  <  j-.  , 

=  J  _  -  2  ’ 

\l-x,  \  <  x  <  1 


.  Then  v 


=  2.  but  v  =  1 


H 


li  -  o 

j  *for  j=l , . . . N 


°°,E,A 
has  | |v|  j 

H 


°°,e  ,A 


However,  any  v  with  L*v 
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We  now  prove  an  imbedding  theorem. 

Lemma  3.2.  Let  v€H2  =  (H2  J  I  •  I  I  ,  )  ,  then  v  CL  (DflH^l)  with 

00  p 

H  . 
p  »e ,  A 


(3.13) 


|v|lLJI)  i  4i  Ih2 

p,E  ,A 


and 


(3.1A)  I  )  v 


'  'l.P.I  - 


1  -1 

P  — 

(1-1) 

-  •-  P 

.  .  2 (1-e  ) 

1 

<  5+  i 

1  .1 

( 

pP(l-e  r) 

J 

p,e  (A 


for  1  <  p  < 


Furthermore, if  we  assume  0  <  e  <  e  ,  then 

o 


(-  -  1) 

(3.15)_  | |v| l1>pjl  1  cc  P  I ! v I  I  2 

Hp  ,e  ,A 


*  ce 


H 


p,e,A 


,  1  ^  1 

1<P<£*>  —  +  — 

-  -  p  q 


=  1  , 


where  c  is  independent  of  r  ,  p  ,  and  A  . 


Proof:  We  use  (3.6)  for  v(x)  and  write 


v(x) 
w1  (x) 


w^(x)  +  w^(x)  +  w^(x)  •  where 
x  _  (x-t) 

f  ( (L*v) (t ) ) ( 1-e  F  )dt  , 


O 
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For  p  *  ®»  we  have 


2l 


•1 


w 


3’  '  L  (1)  -  _  I 

1-e  r 


1.  . 

“V- 


is  the  limiting  estimate  for  p  * 
Next , 


x-x . 


w^,(x) 


l 


x^x 


e  C  J(v’(x.)) 


N-l 

y  zi(x> 


i=i 


where 


0  . 

,  ,  ,  (x-x  ) 

Zj (x)  =  <  i 


if  x^  >  x  , 


f  J(v'(Xi))  ,  if  x.  1  x  . 


Therefore,  we  have 


O.U)  IKIIL(1)  i  TlU,llL  i  \  (i)p|J(V(Xi,)|  i 

D  t—  1  P 


I  -i  I  i  -i 


<  ep  y  e  | J(v*  (x^))  |p?  i  I  M  I  2 


i =1 


P,e,A 


1  <  p  1 


Finally,  wj (x)  = 


x  -  r2L-I-> 

f  1  f 

j  (L*v)(t)  ^  e  dt  = 


(L*v)  (t )  -J-  -)dt 


where 


-x 


<Mx) 


if 

x  >  0 

if 

x  <  0 

which 

* 
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Since  L*v€L  (I)  and  -  <)>(—)  €  L  (I)  ,  Young’s  inequality  implies  that 

d  f;  c  1 


w'  €L  (I)  and 

1  pv 


(3.18) 


K'!l  (i)  -  ML*vilL  (d11^  (i) 

p  p  p 


lL*vl  I,  (n  1  I  M  I  2 

p  "  . 

p,r,A 


,  1  <  p  < 


Combining  (3.1b),  (3.17),  and  (3.18)  with  (3.13)  yields  (3.14). 

Lemma  3.1,  equation  (3.5)  and  a  simple  verification  of  (2.3)  imply  that 
B&(u,v)  is  (l,l)-proper  on  H°  &  *  Hq  f  /\  f°r  ^  <  p  <  ,  where  ^  ^  =  1 

and  ( 1 , ^-proper  on  H*i  A  (  A  •  Also,  for  1  <  p  *  <»•  ,  the  spaces 

11°  .  and  .  are  reflexive,  and  so  in  this  case  we  may  apply  Theorem  2.1 

p, A  P,f  , A 

which  leads  to  the  existence  of  either  the  solution  u  ,  for  the  problem 


BA(u,v)  «  f(v),  Vv€Hp>r>A  , 


or  the  solution  v  of  the  adjoint  problem 

Ba(u,v)  =  f(u),  yu€H“(A  . 

Existence  for  the  cases  p  =  1  or  p  =  <»  does  not  follow  from  Theorem  2.1 
directly,  but  has  to  be  dealt  with  separately. 

—^Young's  inequality  states  that  for  1  5  p  <  ■»  ,  if  s€Lp(ERn),  g€L*(IRn) 
then  h  ■  s*g  exists  a.e.  and  belongs  to  Lp(IRn)  and 

l|h||L  1  I|b||l  ||g||L  . 

p  p  1 

Of  course  to  apply  this  result  to  our  case,  we  should  extend  L*v  by  zero  to 


all  of  IR  . 
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Figure  1:  The  basis  functions  A  and  A  • 

E  a 

Figure  1  shows  the  shapes  of  the  basis  functions  for  the  spaces  S2  A  and  ^ 

Remark:  (4.1)  is  the  space  of  L*-splines  and  (4.2)  was  introduced  In  [5], 
110H12],  [15]— [17].  For  other  various  analyses  of  the  use  of  L*-splines,  we 
refer  the  reader  e.g.  to  [7]— [9] ,  (131. 

We  will  now  study  the  inf-sup  condition  of  the  bilinear  form  B^(u,v) 

°n  S1,A  X  S2,A’  Wh6n  S2,A  =  S2,A  defined  by  (A'2)- 

Theorem  4.1.  Let  S,  ,CH°  ,  s!:  „CH2  „  ,  be  defined  bv  (4.1)  and 

l, A  p,A  2,A  q,e,A 

(4.2)  respectively,  with  ^  ^  =  1  and  1  H  P  i  00  •  Then  the  bilinear  form 

defined  by  (3.4)  is  (C^ ,  C^-proper  with  =  1  ,  and  =  2 

g 

Proof.  Given  v€S„  ,,  select  u  €S,  .  such  that 
2 ,  A  v  1 ,  A 
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a  a  _  a 

uv(x.)  =  -eP|J(v' (x.)) |Po .  PsgnJ(v'(x  )) 
J  J  J  3 


j  =  1. 


,n-l 


for  1  <  p  £  00  (or  1  q  <  oo)  . 


Then  we  have 


N-l  1  + 


(4.4) 


BA(uv,v)  =  lc  P|J(v'(x,))| 


(1  +  q/p) 


P  . 


-q/p 


j  =  l 


N-l 


=  l  eq|J<v'(x  ))|V 
j  =  l  3  3 


q"1'<'  *  I  Ml", 

H 


q,r.,A 


For  any  piecewise  linear  u  we  have 


| u | Pdx  *  l  |u|Pdx  £  l  4( |u(x ,  .)|P  +  |u(x  )|P)h.  < 

J  j=1J  j-12  J"3  3  3~ 


N-l 


i  I  |u(x  ) |Po  1  <  p  < 
j  =  l  3  3 


and  u 


L<*>^  1  <J  <N—  1  }  u  ^  X  j  ^  1 

Therefore, 


<4-5)  lluvllP0  i  2  I  luv(*i> |PP 4  i  2  l  eq|J(v’(x  )) |V 
»p,A  l*1  J  j  J 


N-l 

p~  .  o  v  =  2 1 1 v |  j q , 

H 


q,r  ,A 


i  i  p  n  '■  *  and  1 1  u 


V  H° 
^.A 


Uv ' ' L  (1)  =  1  ’ 


Finallv,  (4.4)  and  inplv 

I  I..I  I1! 


B.  (u  ,v)  H  . 

T^-n - orr-n^-  i  2'l/pIMIV/p  ■  2‘1/pl 


JKJT 

11  vM  o 

H  . 
p,A 


.  1  <  P 


V  H° 
p,A 


H 


q,f.  ,A 


H 


q  ,r  , A 


19 


If  p  =  1  ,  q  =  00  take  u  (x.)  =  d,  ,  where  d.  j  =  is  defined  by 

B  (U  v)  V,  J  J  ’ 

/  '1  t  o  \  mL _  lAV  -Lflll  i  j  i  t  .  .  <  i 


(3.12).  Then 


rn —  "2 '  * v '  *  2 

V  H°  A  H~,£,A 

p,& 


,  which  corresponds  to  the  above 


estimate  when  p  =  1.  Hence,  Theorem  4.1  follows  with  =  2  P, 1  <  p  < 
Obviously  =  1,  and  condition  (2.3)'  is  also  easy  to  check. 

This  leads  to  the  following  theorem. 


Theorem  4.2.  Let  ^  and  ^  be  the  spaces  defined  by  (4.1)  and 

o  F 

(4 . 2),  respectively  and  1  <  p  <  °°  .  If  u€H  and  we  denote  bv  u/€S.  ,  the 

—  —  p  ,A  hi, A 


function  such  that 


Ba(u^,v)  =  B(u,v)  V v €  A 


(4.6) 


luh-ull  0  i  (1+2P)  inf  | |u-w( j 

11  nu  ^ 


Proof.  Theorem  (4.2)  follows  immediately  from  Theorem  4.1  and  Theorem  2.2. 

E 

The  next  theorem  shows  that  (u-u,)(x.)  *  0  ,  for  i  =  1,...,N-1  . 

h  J 


Theorem  4.3.  Let  u€H  ,  1 


<  p  <  °°>  and  u,  €S  €H°  , 
—  —  n  l ,  A  p,; 


,  with 


(4.7) 


B(uh ,v)  =  B<u,v)  ,  V  v  €  sij  A 


(4.8) 


“h (xi)  ~  d£  =  0  ,  i  =  1 . N-l 


with  u  (u,dj , . . .dN-1) 


Proof.  Let  1  <  p  < 


CO 


such  that 


.  By  Theorem  2.1,  there  exists  v,€H 

i  .  p,c,A 


(4.9) 


®A(u,Vi)  =  di 


V 


u  €H 


o 

q,A 


(4.9)  yields  -tv’^'  -  vj  =  0  ,  on  every  1^  ,  J(v£(x  ))  =  0  ,  j  ^  i  ,  and 

-  i  . 

Because  v^(x)  is  continuous,  and  ev^  +  v|  =  0,  on  (O.x^)  and  (x^,l), 

2  E 

v.  is  the  Green's  function  at  x  =  x.  .  Obviously  v.  €H  .  and  v,  €S„ 
i  i  i®,c,A  i2,A 

This  implies  that  (4.9)  holds  for  1  <  q  <  <*>  .  Thus  (4.7)  and  (4.9)  imply 

Bfu^-u.v^  =  0  *  (u^(x1)-d1)  , 

which  finishes  the  proof. 

Theorem  4.3  is  a  restatement  of  the  well  known  fact  that  when  the  Green's 
functions  at  x  =  x^^ ,  1  <  i  <  N-l  ,  belong  to  the  test  space,  then  the  error 
at  the  nodal  points  is  zero. 

Now  we  turn  to  the  study  of  the  space  S2  ^  in  a  more  general  framework. 
Since  B^(u,v)  depends  only  on  the  derivatives  of  the  test  functions,  we  let 

N 

(4.10)  v'  =  l  aX, 

j-1  J  3 

where  Xj  has  support  in  1^  . 

The  situation  is  clear  from  Fig.  2 
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OjXj 


Figure  2:  The  graph  of  v'  =  £  a.x., 

j  =  l  J  J 

1 


where  v€S~  ,  ,  and  uCS 
2,  A’ 


1  .A 


Obviously  we  have 


v'dx  =  0,  because  v(0)  =  v(l)  =  0 


Let  Xj  .  j  =  have  the  following  properties. 


x . 
J 


(4.11)  i)  |  Xj  (x)dx  =  Y^V  , 


ii) 


Xj-1 


texJ+Xjl^dx=  YpV 


j-1 


iii) 


x-x. 


*J-1 


j 


iv) 


Uxj+Xj  |Pdx  * 


'j-1 
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v)  a)  =  Y5  , 


b)  xj(xj)=v;jJ  • 


For  ii  €  S,  A  with  u(x.)  =  d.  ,  j=-l,...N-l,  and  v€  S  ,  (A. 11),  (3.4),  and 


J  .1 


(3.9)  imply 


(4.12) 


Ba(u,v)  =  -_Iidj[a.Y‘JJh.+a.+1y^  +  11hj+1 


[j+il_.  „  Ml- 


~  jl1£dj Iaj+lV5  “ajY6 


(4.13)  a)  |  |v|  |P  =  I  |a  |PyJJV  +  IcP|a  ^5 

h;  „  A  j-i  J  J-* 


[j+i]  _  .  JJIip  i-p 


P  ,g  ,  A 


b)  J  a.Y^h.  =  0  . 

j  =  l  3  3 

In  the  argument  that  follows  we  take  p  =  2,  h.  =  h  for  i  =  ],...N,  where  h  =  ,  and 

Y^  =  Y^(h,t),  j  ~  1,...,N,  depends  in  general  on  h  and  t  .  Before  elaborating  on 
the  next  theorem,  we  will  prove  some  lemmas. 


Lemma  4.4.  Let  a ^ ,  j  =  be  such  that 


(4.14) 


aj~ViK  *  l’  3  =  1,*"’N'1  ’ 


y  aj  =  °  . 

j-i  3 


(4.15) 


V 
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then 


(4.16) 


(a)  <fr(K,N) 


N 

J 

i  =  i 


2 
a  . 
J 


N  2  ^+1 _ 

— - 7:  +  N  » *  /> 

(l-K)2  (KN-1)(K2-1) 


and 


(b)  <KI,N) 


N 

l 

j=l 


2 
a . 
1 


T 

N  -N 
12 


K  4  1 


Proof.  From  (4.14)  It  follows  that 


a .  - 


1 

l-K 


+  AK 


-.1+1 


j  -  1 , .  .  . 


K  t  1 


Using  (4.15)  we  find  that 


and  (4.16)  follows  easily. 

If  K  =  1,  a.  =  ”■  -  j  which  implies  4.16(b)  . 

Lemma  4.5.  Let  <t>(K,N)  be  defined  by  (4.16),  then 


1 


(4.1?) 


*<K,N) 


N(N-l) 

~  U+!k|)2 


Proof.  First  suppose  K>0  (K^D 


Then 


■fettMMfeMtotifcflftittft 


■ia*i  m.i  .iifnlmi 


J 
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*(K,N) 


N  x _ N(1+KN) _ 

(1-K) 2  L  (l-K)(l+K+...+KN-1) 


=  N  N-1  +  K(2N-4)  +...+  KN"2(N-1) 

(l+K)(l+K+...+KN_1) 

2 

Observing  that  jN  -  j  N-l,  for  j  =  1,...,N-1,  and  since  we  are  assuming 
that  K  >_  0  we  get 


4>(K,N)  ■>  N 


(N-1)(1+K+...KN~2)' 

(1+K) (1+. . .+KN~1). 


> 


N(N-l) 

( l+K) 2 


If  K  1  (4.17)  follows  directly  from  (4.16b).  Next  suppose  that  K  <  0  .  If 
N  is  even,  from  (4.16)  we  obtain 


*(K,N>  > 

~  |k2-i| 


N  >  N(N-l) 
(1-K)2  ~  (1+|K|)2 


If  N  is  odd,  then  we  have 


$(K,N)  =  -  - 


(1-K) ' 


N2(KN-1-KN  2+...+l) 
(KN-1)(K-1) 


N  +  N2(KN~1-KN~2+. . .+1) 
(1-K)2  (K-l)2(KN_1+KN~2+...+l) 


Because  N  is  odd  and  K  <  0,  we  have 


kn-1-kn-2+...+i 


kn-x+. 


.+1  >  0  , 


and  so 


♦(K,N) 


(K-l)‘ 


N 

(1-K)2 


N(N-l) 

(i+|k|)2 


which  completes  the  proof  of  the  lemma. 
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Define  \>  and  p  by 


(4.19) 


v  =  y3  -  y6  •  U  =  Y2  +  f  Y5  ,  and  let  K  =  -  * 


2  h  1 5 


when  v  ^  0  . 

Then  we  have,  by  (4.12), 


N-l 


N-l 


(4.20) 


B  (u,v)  =  -h  }  d  [a  .v-Fa,  .  p]  =  -hv  ][  d.[a.-a.  K],  provided  v  i  0 

n  J  =  1  J  J  j—1  ^  ^  -J  '  ^ 


Furthermore,  using  (4.13)  we  have 


(4.21) 


H 


2,e  , A 


1  hY4  l  aj 

j  =  l  J 


Select  v^,  where  v  is  of  the  form  (4.10),  such  that  it  is  characterized  by 
N 

( a .  ■ . _ ^  of  Lemma  4.4.  Then  (4.20)  and  (4.14)  imply  that 


(4.22)  sup  | B  ^ ( u , vq ) |  = 


M  <1 

U°  - 

”2 ,  A 


N-l  -r-  N_i  o  1/2 

sup  |h] |v|  l  d  <  |h]  I v I (  y  dt) 

Mi  n  51  j=l  J  1=1  J 


u€S 


1,A 


u€5 


2, A 
l.A 


where  we  used  the  fact  that 


(4.23) 


7  N_1  ? 

M u  I!  o  l  h  ).  dx 

”2 , A  1=1  J 


Using  (4.21),  (4.22),  and  Lemmas  (4.4)  and  (4.5),  we  obtain 
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(4.24) 


inf 

Ml  2  <} 

H2  ,t  ,A 
u€S2,A 


sup  |B  (u,v)|  < 


1  U i  |  <1 

I  I  «  Q-- 

H 

o 


[4>(K.N)N_\4]2 


v€S 


2, A 


3  1  v  1  ( 1+ 1 K  1 )  ,  3(|p1  +  M) 
_1  1  -  I! 

2  2  2  2 
(N-l)  (N-l) 


A  slight  modification 
^€^(0,1)),  i  =  1, 

(4.25) 


for  v  =  0  will  yield  the  same"  result.  Assume 
,.,m,  m  >  1  ,  are  linearily  independent,  and 


that 


where  depends  in  general  on  c  and  h.  We  then  have  the  following  lemma. 

Lemma  4.7.  Let  c  =  (c^,...cm)  . 

Assume  j-  £  n  is  sufficiently  small.  Then 


(4.26) 


where  L 


sup 

c 


is  independent  of  e 


(|v|2+M2) 

Y4 

and  h  . 


<  ],  <  oo 


Proof.  First  note  that  (4.25)  and  the  definitions  of  v,  ii ,  and  y^ 

2  2 

(see  (4.19)  and  (4.10)  iv  ),  imply  that  Iv)  ,  I  u  |  ,  and  y  ^  are  all 

quadratic  iorms  in  c  =  (i'j,...cm).  Therefore  we  can  restrict  ourselves  to  the 
m  t 

case  where  =  1  . 

The  definition  of  y ^  ,  (see  (4.10)iii)  and  rescaling  imply  that 
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y 


3 


m 


(y  )+4>i  (y)  ]-y  dy 


m 

<  max(l,^)  l  I c  I 
1  =  1 


m 

Also  (fc  =  x.(Xj)  v  h  JJcJ  .  Thus  v  =  Y3-  ^  Y6  1 

m 

max(l,3  I  I c  j  I  !  I  II  ^  2  I  '  °^ta^n  the  same  approximation  for  p  and 

t  ,2m 

so  |v|‘  +  |p|“  <  2  max(  1 ,9  ~)  \  ||<k  ||  . 

h  i=l 

_  2  2 

Since  it  is  assumed  that  ^  is  hounded,  |v|  +  |p|  is  bounded  as  well, 

independent  of  r  and  li  . 

Because  are  linearly  independent  in  L2(I)  and  iflCL^d) 

,  t  m 

i=l,...,m  ,  there  exists  an  n  >  0  such  that  ^  £  ')  implies  that  +—  d> ^ ^ 

are  also  linearly  independent  in  L2(I)  .  The  only  way  for  the  left  hand  side 
of  (4.26)  to  go  to  infinity  would  be  for  y^  to  go  to  zero.  That  is,  there  must 
exist  a  sequence  cf  such  that  y^Ccg)  "*0  as  t  -+  °°  .  But  |cf|  =  1  and 
thus  {cf}  has  a  subsequence  converging  to  say  c^  .  Since  y ^  depends  con¬ 
tinuously  on  c  ,  Y^(cq)  =  0  ,  contradicting  the  linear  independence  of 
{(J>^  +  ^  d>  ^ ^  ,  and  thus  concluding  the  proof. 

Now  we  can  formulate  the  main  result  of  this  section  as  the  following  theorem. 

Theorem  4.8.  Let  h^  =  h  ,  j  =  1,...,N  ,  and  ^  1  h  be  sufficiently  small. 
Let  S2  an  ^  *  dimensional  space,  such  that  if  v€’S2  then 


v '  I  (  =  a^Xj (x)  with 


x-x 

<  (x)  =  l  ci  ( - ^  ) 

J  1=1  1  nj 


x€lj  ,  where 


JC 
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<*iCi 


are  linearly  independent  in 


inf 

Ml  2  =1 

H2,e,A 
v€S, 


sup  <  KN 

HI  0=i  ’ 

HA 


2,  A 


u€S 


1,A 


L2(I) 

JL 

2 


.  Then 


K  >  0  (independent  of 


c  and  h)  . 


Proof.  The  theorem  follows  immediately  from  (4.24)  and  Lemma  4.7. 

It  is  readily  seen  that  the  space  ^  satisfies  the  assumptions  of 
Theorem  4.8  for  any  a  .  Theorem  4.8  shows  that  the  use  of  test  functions 
recommended  in  [5],  [10]- [12],  [15]- [17],  leads  to  an  inf-sup  condition  converging  to 
0  as  N  *  ,  when  ^  is  held  constant.  Moreover  ,  there  is  no  choice  of  4>^(x)  , 
independent  of  e  and  h  ,  which  would  lead  to  quasi-optimal  results.  Neverthe- 
less^in  the  next  section  we  will  show  that  the  situation  significantly  improves 
when  assumptions  about  the  input  data  are  made. 
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5.  MODIFIED  VARIATIONAL  PRINCIPLE 

In  this  section  we  perform  a  more  detailed  analysis  of  the  subspace  ^  , 

defined  by  (4.3).  In  particular,  we  will  compare  this  space  to 
E 

S»2  ^  ,  when  each  is  being  tested  against  the  subspace  of  piecewise 

linear  trial  functions, 

E  ot 

To  distinguish  between  the  spaces  S„  and  S„  ,  we  will  write  v„  and 

Z ,  A  Z ,  A  b 

E  a 

vq  ,  denoting  a  function  in  ^  and  S2  ^  respectively.  By  straightforward 
computations  we  get 


(5.1) 


where 


(5.2) 


VE(x)>I, 


EJ(X)  * 


vrx 


fj  (x) 


h, 

-J-  _S_ 


(l-e-hJ/£)  ’ 


vj  =  vE(xj) 


Analogously , 


(5.3) 


v;(x)| 


v.-v.  . 

h.  xj(x)  ’ 
j 


where 


(5.4) 


x.(x)  =  1  -  T-^-  (x-x  .)  ,  X 

J  j  j'2  j 


VVi 


and  once  again  v.  *  v  (x.)  . 
6  j  «  j 


Let  P  be  the  orthogonal  projection  operator  from  L2 ( I ^ )  to  linear  functic 
in  Ij  .  By  direct  computation  we  obtain 
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(5.5)  P  E  (x)  =  1  -  ^-(coth  ^  =  1  -  t  aj(x~x  1>  *  *°(x)  * 

j  J  j  j-  2  j 


where 


(5.6) 


h  , 

a?  =  a?(e,h  )  =  coth  o  - 
J  J  j  2e  hj 


a  o 

From  now  on  we  denote  A  =  ^2  A  *  that  is,  we  consider  the  particular 
a  =  a°5  where  a°  is  defined  by  (5.6).  We  remark  that  a°  is  the  same  as  the 
"optimal"  a  described  in  [11],  [12],  [15]-[17]. 

E  a 

We  now  present  some  important  relationships  between  ^  and  S2  A  . 

Lemma  5.1.  Let  v  € S,  and  v  € .  ,  (a  =  a°) ,  be  such  that  v  (x.)  = 
E  z,A  a  l ,A  E  j 

v  (x  )  for  j  *  0, . . . ,N .  Then  we  have 
^  J 


(5.7) 


(vE’z)I  =  (va’z)I 


for  any  z€L2(I)  ,  z  linear  on  1^  for  j  =  1,...,N,  (  z  not  necessarily 


continuous  at  Xj),  and 


(5.8) 


(vE’q)I  *  (Vq)I 


for  any  q€L2(l)  ,  q  constant  on  Ij  for  j  =  1,...,N. 


Proof.  (5.7)  follows  immediately  from  (5.1)— (5.6) .  To  prove  (5.8)  let 
q  be  piecewise  constant  on  A  .  Let  z(x)  *  fXq(t)dt  .  Then  z  is  linear  on 
Ij  and  therefore  by  (5.7)  we  have 
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1  1 
v’zdx  =  v'zdx  . 

J  t  J  a 

0  0 

Integration  by  parts,  recalling  that  z  is  continuous,  and  that  v£(0)  = 

v  (0)  and  v  (1)  =  v  (1),  yields  the  desired  result, 
ot  c.  ot 

Theorem  5.2.  Let  u€S.  .  ,  v_  ,  v  be  as  in  lemma  5.1.  Then 

1,A  E  a 


(5.9)  Vu,vE>  =  BA(u’va)  * 

Furthermore,  let  f  be  a  linear  combination  of  Dirac  "delta  functions" 
centered  at  the  mesh  points  x.  ,  j  =  1,...,N-1,  and  any  piecewise  constant 
function  on  A  .  Let  u,,€S,  and  u  €  S.  be  such  that 


and 


Then 


BA^uE’vE^  *  ^f’vE*  ’  ^vE€S2, 

BA(ua>V  *  (f’V  *  *Va€S2,A  * 


<5* 10)  UE  =  ua 

E  a 

Proof.  For  u€S^  &  and  v€S£  ^  or  v€S^  A  we  have  BA(u,v)  =  (eu'-u.v'). 

Since  u'  -  u  is  linear  on  1^  ,  (5.7)  implies  (5.9). 

The  fact  that  v„(x.)  =  v  (x.)  for  j  =  1,...,N-1,  and  (5.8)  together  with 
t  j  a  i 

(5.9)f  yield  (5.10). 

Theorem  5.2  shows  that  the  global  stiffness  matrices  generated  by  using  ^ 
with  Sh;  and  S™  ,  (n=u  ),  respectively,  are  identical.  However,  we  have  shown  (see 
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Theorem  4. 2),  that  for  any  u€H°  . 

p,A 


(5.11) 


£  c 
A 


inf 
w€S,  , 


with  C  independent  of  A  ,  e  ,  On  the  other  hand,  we  have  shown  (see  Theorem 

4.8)  that  there  exists  u€H°  ,  (resp.  f  €  (H~  .)'),  such  that 

Z  •  A  Z ,  c  ,  A 


(5.12) 


on  a  uniform  partition.  Comparing  (5.11)  and  (5.12)  we  see  that  the  use  of 

Z  ,  A 

will  deteriorate  the  result  for  some  u,  (rasp.  £).  The  question  arises  as  to 
whether  or  not  this  f  has  any  "engineering  meaning",  that  is,  will  this  effect 
be  observed  in  practice. 

The  preceding  theorem  says  that  if  f  is  piecewise  constant,  then 
ug  =  ua  •  This  suggests  that  if  f  is  "nearly"  piecewise  constant,  then  u^ 
may  be  "nearly"  equal  to  u^  .  We  shall  now  show  that  this^is  true;  that  is, 
if  f  (the  input  data)  is  "reasonable"  then  the  term  Ch  can  be  replaced 
by  C  independent  of  e  and  h  . 

Theorem  5.3.  Let  f  €  (H^  e  fi) '  ,  fCH^d^)  j  =  ],...,  N.  Denote  by 

u€H  the  exact  weak  solution  of  the  problem 
p  ,A 

-Eu"  +  u’  =  f  , 


u(0)  ■  u(l)  *  0 
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that  is. 


(5.13) 


B  (u,v)  =  (f,v),  Vv€Hq>f^ 


Let  u  .  be  such  that 

a  i » & 


(5.14) 


B.(u  , v)  =  (f,v)  ,  VvCS^  A  ,  a  =  a 
A  a 


Then  for  1  1  P  1  ”  > 


(5.15) 


lu-u  ||  <  (1+21/P)  inf  I  I u-w 1 |  + 

1  -°  "  w€S.  .  H° 


H 


p.A 


i,A 


+  C  max  [tu  1 1  f  j  |  ,  1  » 

j-l . NU 


with  C  independent  of  A  ,  c  ,  p  ,  and  f 
We  first  prove  a  lemma. 

Lemma  5.4. 


|v'-v^l t,  ,T  x  i  cl Iv^l lL  (!)  *  1  -  q  - 

q  j 


^“VE»'Lq(lj, 


with  C  independent  of  q  ,  A,  and  t. 

Proof.  Obviously  it  is  sufficient  to  show  that 


(5.16) 


KI,Lq(Ij)  -Cl,VEilLq(Iq) 


♦ 


Since  v'  is  the  projection  of  vi  onto  piecewise  linear  functions,  we  have 
a  t 


(5.17) 


and  hence 


(vE“V«’Va)IJ  =  ° 


(v*)2dx  = 
n 


I . 
3 


vEVadx  -  * lVE»  lLq(Ij) ! ,V«!  Lp(Ij)  • 


Therefore , 


(5.18) 


[I 

Vi - , , 

|  k;iPd*jP 


|v' |qdx  ,  1  1  q 


(Where,  if  q  =  1,  we  replace  the  denominator  by  I  I i I L  (I  ) 

00  j 


.) 


Returning  to  (5.1)  and  (5.3),  we  see  that  instead  of  and 

can  dea  l  with  Xj  (x)  and  Ej(x)»  resPective3y •  Since 


V  we 


and 


Xj(x)  1  1  +  3aj  ’ 


Xj  Cx)  >  1  , 


*€I3 


X,_!  1  x  1  x  ! 


\J-1 


j- 


we  get 


2  q  h1  q 

X,(x)dx)q  >  (-^-)  , 


'j 


1 


mmmMMBIh 


f 
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and 


(I  l*j|Pdy) 

TJ 


a. 

P(l+3aJ)~q2_q  =  hj(2(l+3aj))"q 


However  , 


|  I Xj  |qdx  <  hj(1+3aj>q  , 


which  implies 


(5.20) 


q 

p  X2(x)dxj 


1  E|qdx  >  - J - -g;  1  h  .  ( (l+3a.)  2)' 

Ij  £  |x(x)(PdxP 


>  2  q(l+3aj)'2q|  | Xj |qdx 


Th  refore,  by  (5.18)  -  (5.20)  we  have 


VV  1 


2(H3aj)2||v'||^(Ij) 


(5.21) 
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Furthermore,  it  is  easy  to  see  from  (5.6)  that  0  <  a  <  1  ,  and  thus 

because  of  (5.21),  (5.16)  holds  for  l^q<°°.  If  q  =  ®  we  simply  note  that 

I  IE, I  I.  .  <  1  +  3a.  and  hence  (5.16)  also  holds  in  this  case.  Now  we 

j  L°°(Ij )  -  J 

will  prove  Theorem  A. 3. 


Proof.  Let  v  (x.)  =  v  (x.)  for  j  =  0,...N  . 

E  j  a  j  J  ’ 


Since 


(f’Va}  =  WV  =  VUa’V  * 


and 


(f,vE)  =  Ba(uf,ve) 


f 


we  get 


B<uE~ua»vE)  =  <f»vE”v0)  • 


Hence,  applying  Theorem  A.l  yields 


!VV 


p,A 


<  C  sup.- 

Vs* 

E  2,  A 


(f’VE"v«) 


H 


q  ,e ,  A 


By  Theorems  4.2,  4.1,  2.2,  and  the  triangle  inequality 
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(5.22) 


u-u  <  u-u_  +  u„-u  < 

a1  1  o  E1  1  o  1  1  E  a1 '„o  — 

H  .  H  H 

p,A  p,&  p , A 


1/p 

<  (1+2  )  Inf 

w€S 

A  » 


I  u-w  M  +  C  sup 

-O 


(f ,V  -V  ) 

E  a 


H'  .  vc€S^  11  E 1  '2 
p,A  E  2, A  H 

q,f 


Therefore,  we  have  to  deal  with  the  term 


sup 


(f*vE"V 


<eE  E 1 1  2 
v  6S„  .  H  . 

E  2, A  q.c.A 


Denote  by  £  the  function  which  is  constant  on  I j  ,  which  satisfies. 


(f-?)dx 


I  . 
3 


If  we  let  F . (x)  * 

3  J 


(f-?)(t)dt  ,  then  Fj  (Xj  i)  =  Fj(x^)  =  0,  and  since 


3-1 


f  €H  (I . )  ,  we  have 

CO  J 


(5.23) 


I Fj (x> I  i  hj I I li.co.i 


j 


Further,  we  denote  by  F(x)  the  function  such  that  F(x)  |  =  F.(x)  .  By 

J  J 

Lemma  5.1  we  get 


=  0  ’ 


(f,vp-vrt)  «  (f-f  .Vp-vJ  =  (F* » v£*v^) 


and  hence 
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Integration  by  parts,  (5.23),  Lemma  5.4  and  Lemma  3.2  yield 


i<f.vv>i  i  .n<vvE-v;>i .i ; iiv'-v;ii  ; 

i_i  j  j-i  j  i  j 


<  c 


E  11  2 


l.e  ,A 


It  can  be  easily  verified  from  (3.9)  and  (4.2)  that 


1 

|vEll  2  12’llvjll  2  ,  liq; 

q.c.A 


l  has,  we  have 

(5.24) 


SUE 


.  1.2 


<  C  max  h .  f 


v  CS^'  „2 

®  2»A  ^q,e,A 


hJl lf 111  t  1  • 


with  C  independent  of  A,  q,  e  and  f.  (5.24)  together  with  (5.22)  proves  the 


theorem. 
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6.  NUMERICAL  CONSIDERATIONS 

Theorem  5.3  shows  that  for  general  meshes  and  parabol ically  upwinded  test 
functions,  the  error  consists  of  two  parts.  In  tnis  section  we  comment  on  the 
interplay  between  them. 

for  the  error  estimate  (5.15)  to  be  "nearly"quasi-optimal ,  the  second  term 

(6.1)  C  max  [h2 | | f | |  1  ] 

j=l,...N  J  H»(lj) 

should  be  smaller  than,  or  of  the  same  magnitude  as,  the  first  (quas i-opt imal) 
term 

(6.2)  (l+21/p)  inf 

w€  S 

l,u 

It  Is  well  known  from  singular  perturbation  theory  that  the  exact  solution 
to  (1.1)  can  be  written  as  a  combination  of  a  smootn  term  and  a  boundary  layer  term. 
The  boundary  layer  term  is  of  the  form  e^X  ^  C  ;  thus,  it  is  sufficient  to  analyze 


inf 
w€S 

A  simple  analysis  of  this  term  on  a  uniform  mesh  with  stepsize  h  implies  that 
for  c^sh 


(x-l)/r  m 
e  -w||Ho 

p,a 
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to  (6.2). 

If  t>>h,  standard 

2 

h  ,  the  same  magnitude 
Finally,  we  remark 
In  a  forthcoming  paper, 
will  be  studied. 


approximation  theory  says  that  (6.2)  will  have  magnitude 
as  (6.1). 

that  the  constant  C  in  (6.1)  can  be  computed  exactly, 
an  adaptive  approach  with  a-posteriori  error  estimates 


41 


REFERENCES 


BabuSka,  I.,  Aziz,  A.K.,  "Survey  Lectures  on  the  Mathematical  Foundation 
of  the  Finite  Element  Method",  in  The  Mathematical  Foundations  of  the 
Finite  Element  Method  with  Applications  to  Partial  Differential  Equations , 

A.K.  Aziz,  ed. ,  pp. 5-359.  New  York:  Academic  Press,  1973. 

Babuska,  I.,  Osborn,  J.,  "Analysis  of  Finite  Element  Methods  for  Second 
Order  Boundary  Value  Problems  using  Mesh  Dependent  Norms",  Numer.  Math. 

34 ,  41-62  (1980). 

Barrett,  J.W.,  Morton,  K.W.,  "Optimal  Finite  Element  Solutions  to  Diffusion- 
Convection  Problems  in  One  Dimension",  Int.  J.  Num.  Meth.  Engng. 

(1980) 

Barrett ,  J.W.,  Morton,  K.W. ,  "Optimal  Petrov-Oalerkin  Methods  Through 
Approximate  Symmetrization" ,  U.  of  Reading,  Num.  Anal.  Report  4/80. 

Christie,  J.,  Mitchell,  A.R.,  "Upwinding  of  High  Order  Calerkin  Methods  in 

Conduction-Convection  Problems",  Int.  J.  Num.  Meths.  in  F.ngng.  1_2 ,  1  764- 

1771  (1978).  '  \ 

Christie,  J.,  Griffith,  D.F.,  Mitchell,  A.R.,  Zienkiewicz,  O.C.,  "Finite 
Element  Methods  for  Second  Order  Differential  Equations  with  Significant 
First  Derivatives",  Int.  .1.  Num.  Meths.  in  Engng.,  10,  1389-1396  (1976). 

DeCroen,  P.P.N.,  "A  Finite  Element  Method  with  Large  Mesh-width  for  Stiff 
Two-point  Boundary  Value  Problems",  Preprint,  Eidenhoven  University, 

Department  of  Mathematics,  Eidenhoven,  The  Netherlands  (1978). 

DeGroen,  P.P.N.,  Hemker,  P.W.,  "Error  Bounds  for  Exponentially  Fitted 
Calerkin  Methods  Applied  to  Stiff  Two-point  Boundary  Value  Problems", 

Numerical  Analysis  of  Singular  Perturbation  Problems  (eds.  P.W.  Hemker 
and  J.J.H.  Miller),  Academic  Press  (1979). 

Diaz-Munio,  R.F.,  Wellford,  C. ,  Jr.,  "A  Finite  Element  Singular  Perturba¬ 
tion  Technique  for  Convection-Diffusion  Problems.  Part  I:  The  One- 
Dimensional  Case",  to  appear  in  ASME  J.  of  Applied  Mechanics. 

Griffiths,  D.F.,  Lorenz,  J.,  "An  Analysis  of  the  Pet rov-Galerkin  Finite 
Element  Method",  Comp.  Meth.  Appl.  Mech.  Engng.,  1 4 ,  39-64  (1978). 

Heinrich,  J.C.,  Zienkiewicz,  O.C.,  "The  Finite  Element  Method  and  Upwinding 
Techniques  in  the  Numerical  Solution  of  Convection  Dominated  Flow  Problems", 

Finite  Element  Methods  for  Convection  Dominated  Flows,  ed.  by  J.R.  Hughes, 
pp. 105-135  (1980). 

Heinrich,  J.C.,  Huyakorn,  P.S.,  Zienkiewicz,  O.C.,  Mitchell,  A.R.,  "An 
"Upwind"  Finite  Element  Scheme  for  Two-Dimensional  Convective  Transport 
Equation",  Int.  J.  Num.  Meth.  Engng . ,  11,  131-143  (1977). 


13.  Hemker,  P.W.,  A  Numerical  Study  of  Stiff  Two-Point  Boundary  Problems, 
Mathematisch  Centrum,  Amsterdam  (1977). 


42 


14.  Hughes,  T.J.R.,  Liu,  W.K. ,  Brooks,  A.,  "Finite  Element  Analysis  of 
Incompressible  Viscous  Flows  by  the  Penalty  Function  Formulation",  .J. 

Comp.  Physics  30,  1-60  (1979). 

15.  Mitchell,  A.R.,  Christie,  J.,  "Finite  Difference/Finite  Element  Methods 
at  the  Parabolic-Hyperbolic  Interface",  Proceedings  Conf .  Numerical 
Analysis  of  Singular  Perturbation  Problems,  Catholic  University,  Nijmegen, 
Holland,  May  1978.  Academic  Press,  London  (1979). 

16.  Mitchell,  A.R. ,  Griffiths,  D.F.,  "Upwinding  by  Petrov-Galerkin  Methods 

in  Convection-Diffusion  Problems",  Report  NA/36,  Department  of  Mathematics, 
University  of  Dundee  (1979). 

17.  Zienkiewicz,  O.C.,  Heinrich,  J.C.,  The  Finite  Element  Methods  and  Con¬ 
vection  Problems  in  Fluid  Mechanics,  Finite  Elements  in  Fluids,  Vol.  3 
(Gallagher,  Zienkiewicz,  Oden,  Morandi-Cecchi  and  Taylor,  eds.),  Wiley, 
London  (1978). 


